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SIGNIFICANCE  AND  EXPLANATION 


Mathenatical  optimization  techniques  are  now  widely  used  in  business,  govern- 
ment service,  and  the  armed  forces.  A typical  problem  involves  finding  the  cheapest 
method  of  providing  certain  services,  using  various  available  sources  of  supply, 
under  given  constraints,  e.g.  on  the  level  of  production  of  individual  sources. 

Problems  where  we  know  exactly  the  sources  of  supply  and  their  limitations,  are 
called  deterministic . Such  problems  have  led  to  intensive  development  of  mathemati- 
cal optimization  techniques  such  as  linear  and  nonlinear  programming. 

In  most  situations  there  is  some  uncertainty  as  regards  the  sources  of  supply 
and  their  limitations.  Such  problems  are  called  stochastic.  They  have  received 
much  less  attention  than  deterministic  problems,  partly  because  they  are  much  more 
difficult. 

The  type  of  problem  considered  in  this  paper  will  be  described  in  terms  of 
electrical  power  systems;  it  should  be  clear  that  the  model  covers  very  many  other 
situations.  Suppose  that  the  power  systems  are  interconnected  so  as  to  form  a net- 
work. If  a system  cannot  serve  the  demand  of  its  own  area  because  of  (random)  excess  i 
demand  and/or  (random)  deficiency  in  the  electricity  generation,  then  the  other  sys- 
tems help  to  meet  the  total  demand.  How  can  we  design  the  cheapest  possible  system 
so  as  to  reach  a prescribed  (high)  reliability  level? 

The  problems  considered  in  this  paper  are  based  on  the  following  underlying 
deterministic  model: 

Minimize  (c'u  + d'v) 
subject  to  Au  ^ b,  Tu  + Mv  ^ C 

When  the  ^ is  assumed  to  be  random,  the  model  is  formulated  as  a two-stage  problem 
in  which  we  first  decide  on  u,  then  choose  the  random  variable  5,  then  find  v. 

The  variable  u represents  supply,  5 represents  random  demand  plus  deficiency  in 
various  areas,  and  v represents  dispatch  power. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive  summary 
lies  with  MRC,  and  not  with  the  author  of  this  report. 


PLANNING  IN  INTERCONNECTED  POVIF.R  SYSTEMS: 


AN  EXAMPLE  OF  TWO-STAGE  PROGRAMMING  UNDER  UNCERTAINTY 
Andras  Prekopa 

1 . Two-Stage  Proqrainininq  under  Uncertainty 

Vectors  and  matrices  will  be  denoted  by  small  resp.  capital  latin  letters.  As  in  many 
papers  and  textbooks  on  probability,  remdom  vectors  will  be  denoted  by  qreek  letters.  The 
prime  means  transpose. 

The  stochastic  programtiing  model:  two  stage  programming  under  uncertainty  has  been 
introduced  by  Dantzig  and  Madansky  [3] . Before  formulating  this,  we  remark  that  stochastic 
programming  models  are  usually  formulated  in  such  a way  that  first  we  start  from  a deter- 
ministic problem  that  we  call  the  underlying  deterministic  problem,  then  observe  that  some 
of  the  parameters  appearing  in  the  model  are  random  and  finally  formulate  a decision  princi- 
ple taking  into  account  the  underlying  problem  and  the  probability  distribution  of  the  ran- 
dom variables  involved. 

In  case  of  the  model  formulated  by  Dantzig  and  Madansky,  the  underlying  deterministic 
problem  is  the  following 

minimize  (c'u  -t  d'v) 
subject  to 

(1.1)  Au  > b , 

« 

Tu  + Mv  > C . 

Having  observed  that  5 is  random,  the  model  of  two-stage  programming  under  uncertainty  is 
formulated  as  the  following  non-linear  programming  problem 

minimize  {c’u  + Ely(u,t)]l 
subject  to 

(1.2)  Au  > b , 

B 

u € K , 
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where  K is  the  set  of  those  u vectors  which  have  the  property  that  for  every  t in  the 
set  of  its  possible  values  there  exists  at  least  one  v (dependinq  on  u and  t)  such  that 
Tu  + Mv  > C and  w(u,C)  is  the  optimum  value  of  the  following  linear  programming  problem 
(where  u and  C are  fixed  for  a moment,  u f K,  Au  > b) : 

minimize  d'v 
subject  to 

(1.3)  Tu  + Mv  > f.  . 

This  is  called  the  second  stage  problem  cind  (1.2)  the  first  stage  problem.  In  the  practical 
application  of  this  model  first  we  decide  on  u on  the  basis  of  Problem  (1.2),  then  t is 

realized  and  finally  we  decide  on  v on  the  basis  of  Problem  (1.3).  We  solve  the  first 

stage  problem  only  once  but  then  the  "experiment"  producing  the  random  varieible  C is 
"performed"  many  times  independently  of  each  other.  Thus  we  solve  the  second  stage  problem 
many  times  and  E(U(u,^))  is  interpreted  as  the  arithmetic  mean  of  the  obtained  optima. 

If  C has  a discrete  distribution,  then  we  can  obtain  at  once  the  first  stage  and  the  second 
stage  optimal  solutions  (3). 

We  recall  a few  facts  concerning  this  model.  We  also  give  the  main  ideas  of  the  corre- 
sponding proofs.  For  more  detailed  explanation  and  references  the  reader  may  consult  the 
boo)(  of  Kail  [6]  . 

K is  a convex  polyhedron  (141  . In  fact  letting  C vary  in  the  entire  space  from 

where  j t2U(es  its  (vector)  values,  the  set  of  all  those  u,5  vectors  for  which  a v can 

be  found  satisfying  Au  + Bv  ^ 5 , forms  a convex  polyhedron.  Assume  that  it  is  given  by 
the  linear  inequalities 

G £ < Hu  . 

Now,  if  h is  the  vector  consisting  of  the  components 

sup  GC 

where  the  supremum  is  talten  with  respect  to  all  realizations  of  ^ (a  zero-probability  set 
can  bo  neglected,  of  course),  then  clearly 

K - (u|hu  j h)  ; 

this  proves  the  assertion. 
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For  ovory  fixed  ^ the  seeond  staqe  optlmam  value  ij(u,f.)  is  a convex  function  of 
the  vai lable  u.  This  follows  from  the  fact  that  the  second  constraint  in  Problem  (1.3) 


can  be  written  as  Mv  ^ f,  - Tu  and  that  the  optimum  value  of  a linear  programming  problem 
is  a convex  function  of  the  right  hand  side  vector  provided  we  have  ^ul  objective  function 
to  be  minimized.  If  the  set  of  [xrssible  values  of  C is  convex, then  p(u,C)  is  convex 
both  in  u and 

If  E(f.)  exists  and  ki(u,C)  is  finite  for  every  feasible  u and  every  t in  its  set 
of  possiiile  values,  then  E(ii(u,C))  exists  and  is  a convex  function  of  u.  In  fact  the 
dual  of  Problem  (1.3)  contains  the  com[x>nents  of  C as  coefficients  of  variables  in  the 
objective  function.  If  such  ^ul  objective  function  is  finite  for  every  C and  E(C)  exists 
then  - as  it  is  easy  to  see  - the  expectation  of  the  objective  function  also  exists.  The 
convexity  of  E(u(u,t))  follows  from  the  convexity  of  li(u,C)  for  every  fixed  f.. 

In  (‘31  we  gave  a variant  for  the  model:  two  staqe  programming  under  uncertainty  by 
dropping  the  condition  that  the  second  stage  problem  must  have  feasible  solution  for  every 
f..  We  require  the  solvability  of  the  second  stage  problem  only  by  a certain  probability. 

The  formulation  of  this  model  was  motivated  by  the  fact  that  in  many  applications  (e.g.  in 
problems  of  engineering  design)  the  working  of  the  system  cannot  be  ensured  by  probability  1 
but  only  by  a somewhat  lowi'r  probability.  In  other  words,  in  many  cases  Problem  (1.3)  will 
not  have  feasible  solution  for  every  [xissible  C-  ff  e.g.  C has  a nondegenerated  multi- 
variate normal  distribution,  then  the  set  of  possible  values  of  J is  the  entire  space  and 
the  requirement  of  the  solvability  of  the  second  stage  problem  for  every  t limits  the 
structure  of  the  matrix  M a great  deal. 

The  sot  of  all  u,^  vectors  for  which  there  exist  v satisfying  the  constraints  of 
Problem  (l.l),  can  be  given  by  a system  of  linear  ineijual  it  ies  with  variables  u and 
Let  us  write  it  (as  before)  in  the  following  form 

(1.4)  G C < Hu  . 

Then  our  new  first  stage  problem  is  the  following 
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minimize  {c'u  + E(ji(u,5)l  I 


1 


subject  to 

Au  > b , 

P{G  C < Hu)  > p , 

where  p is  a prescribed  (high)  probability  and  (u,t)  is  the  random  optimum  value  of  the 
new  second  stage  problem: 

minimize  (d'v  + t'z) 

(1-5)  subject  to 

Tu  + Mv  + z > £ 

Here  u is  fixed  and  such  that  it  satisfies  the  constraint  Au  > b.  The  term  d'v  in  the 
objective  function  expresses  the  cost  of  the  second  stage  activity  and  t’z  is  the  cost  of 
infeasibility  of  the  original  second  stage  problem  (1.3).  We  assume  that  t is  chosen  in 
such  a way  that  all  of  its  conpcments  are  positive  thus  maicing  z = 0 whenever  to  the  given 
u and  £ there  exists  v such  that  Tu  + Mv  > £. 

In  our  problem  of  power  systems  planning  we  shall  use  separeible  nonlinear  function 
I instead  of  the  linear  function  t'z.  It  is  not  difficult  to  show  that  if  the 

functions  non-decreasing  and  convex,  then  lJ(u,£)  is  still  convex  in  u for 

every  fixed  £. 

The  model  of  two-stage  programming  under  uncertainty  with  probabilistic  constraint  on 
the  solvability  of  the  second  stage  problem  will  be  shown  in  this  paper  an  important, 
practical  model  construction.  . 


2.  The  Problem  of  Planning  in  Interconnected  Power  Systems 


Finding  the  reliability  of  interconnected  power  systems,  or  in  other  terms,  finding 
the  Loss  of  Load  Probability  (LOLP)  is  an  important  subject  in  power  systems  engineering 
and  has  an  extensive  literature  [1] . Planning  under  reliability  constraint  seems  to  be  a 
recent  idea.  One  approach  has  been  presented  by  Scherer  and  Joe  [12] . They  assume  that 
each  plant  has  only  one  generating  unit  that  can  be  up  or  down.  Using  mixed  integer  pro- 
greunming  they  find  the  sizes  of  the  plants  subject  to  a probabilistic  constraint.  The  inter 
connection  is  disregarded.  A discussion  of  the  previous  planning  models  is  also  presented. 
Our  approach  will  be  different,  the  interconnection  of  the  power  systems  plays  an  important 
role  here.  Below  we  list  our  assumptions  and  notations. 

Assumptions 

1.  The  power  systems  are  connected  in  such  a way  that  (mathematically)  they  form  a 
networic  (see  the  definition  in  Section  3)  . 

2.  The  system  must  have  a realiability  level  1 - LOLP  > p where  p is  a probability 
(near  unity  in  practice)  prescribed  by  ourselves. 

3.  We  may  increase  the  power  generating  capacities  of  the  existing  systems  and  the 
tie  line  capacities  among  them  for  cost. 

4.  No  distinction  will  be  made  Isetween  capacity  and  reserve  capacity  or  demand  and 
deficiency.  We  shall  simply  use  the  terms  "capacity"  and  "demand". 

5.  Each  power  system  can  help  others  to  the  extent  of  its  surplus  power  generation 
capacity  and  the  tie  line  capacities.  The  tie  line  between  systems  a and  b is  not 
necessarily  the  same  as  that  between  b and  a,  further  these  two  may  have  different 
capacities. 

6.  Power  dispatching  and  outage  cost  money.  The  first  one  will  be  assumed  to  be 
linear  while  the  seond  one  nonlinear  of  a fixed  charge  type. 


Notations 


N number  of  [lower  systems  in  the  pool; 

generating  capacity  of  the  jth  system,  to  be  determined; 

y^l^  capacity  of  the  tie  line  between  system  h and  system  k,  to  be  determined; 

i^  bus  power  at  the  jth  system;  it  can  take  on  positive,  zero  and  negative  values 

depending  whether  plant  j helps  other(s),  serves  only  its  own  area  or  receives 
help  from  others; 

i^j^  power  dispatched  from  system  h to  system  k;  it  is  assumed  that  ij^j^  = - 

|i..  if  i..  > 0 , 

hk  hk  = 


r if  ij,,  <0  , 


0 if  ij,^  > 0 , 


-i.  , if  i.  . < 0 ; 

. hk  hk 


p prescribed  lower  level  for  1-LOLP  ; 

random  demand  for  electric  power  (including  deficiency)  at  the  jth  system; 

Cj (x)  cost  function  of  the  generating  capacity  at  the  jth  system; 

Chk<Y)  cost  function  of  the  interconnection  capacity  between  system  h and  system  k; 

d,  , cost  of  dispatching  one  unit  of  power  from  system  h to  system  k (dispatching 
hk 

cost  is  assumed  to  be  a linear  function) ; 

t^(z)  cost  of  outage  of  magnitude  z at  system  j when  no  other  system  is  able  to 
assist  to  meet  total  demand; 
prescribed  lower  resp.  upper  bound  for  ; 

prescribed  lower  resp.  upper  bound  for  y . 

nK  hK  nK 

Let  A(x,y,C)  denote  the  event  that  for  fixed  x and  y the  total  demand  in  the 
•"ool  can  be  met  by  a suitable  power  dispatch.  The  event  A(x,y,5)  wil]  be  given  in  terms 


of  linear  inequalities  in  the  next  section. 


1 


We  define  further  a dispatch  function  d(x,y,5)  that  is  the  cost  of  the  minimum  cost 

dispatch  necessary  to  meet  the  total  demand  in  the  pool  provided  it  is  possible  i.e. 

ftlxty/f.)  occurs.  In  general,  no  formula  can  be  given  for  d(x,y,5):  for  every  x.y,^ 

the  minimum  cost  dispatch  (flow)  can  be  determined  algorithmically. 

If  A(x,y,5)  does  not  occur,  then  instead  of  the  area  demands  - x^ , j - 1,..,N 

we  use  the  increased  demands  5^  - x^  + , j = 1,...,M,  add  the  sum  of  the 

t (z.),  j = 1,...,N  values  to  the  cost  of  the  flow  and  treat  z,,...,z.,  as  nonnegative 
j j IN 

variables.  Both  from  the  theoretical  and  practical  point  of  view  it  is  desirable  to  chose 
the  outage  costs  in  such  a way  that  t^Cz)  is  monotonically  increasing  in  (0,“)  further 
there  is  a fixed  charge  at  z = 0,  i.e. 


t^(z)  = 


0 if  z = 0 


T.  ^ 0 if  z > 0 
1 


Clearly  we  can  choose  T, ,...,T„  so  large  that  if  A(x,y,£)  occurs,  we  obtain  automati- 

1 N 

cally  z = ...  = z =0  in  the  optimal  solution  of  the  just  described  extension  of  the 
1 N 

flow  problem.  Let  u = u(x,y,C)  denote  its  optimum  value. 

Now  we  formulate  our  planning  model  in  the  following  way 


minimize 


I c.(x,)  + I + E(u) 


1 3 


(2.1) 


j - ^ h,)c 

subject  to 

P(A(x,y,5))  > P , 


X«>  < X.  ^ 

1 =3^3 


(u) 


Y*'’  < V < Y 

h)c  = ^h)i  = h)c 


all  j 

all  h,)i  . 
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3.  Application  of  the  Gale-Hof fman  Theorem  for  the  Determindtion  of  h(x,v,i'.)  in  Terir.s 
of  Linear  Inequalities. 

The  feasibility  theorems  formulated  by  Gale  resp.  Hoffman  concerninq  network  flows 
resp.  circulations  are  essentially  the  same.  The  presentation  given  by  Gale  suits  very  well 
to  our  problem.  Thus  we  follow  his  way  of  thinkina. 

A network  lN,y)  consists  of  a finite  set  of  nodes  N and  a capacity  function  y 
defined  on  the  product  set  N » N,  where  y(a,b)  can  take  on  nonnegative  values  and  + ■». 

A flow  i on  [N,yJ  is  defined  as  a function  on  N « N satisfying  the  relations 

i(a,b)  + i(b,a)  = 0 , 

for  all  a,b  ( N 


i(a,b)  < y(a,b) 


A network  demand  q is  a real  valued  function  on  N.  .A  negative  q will  be  inter- 
preted as  a supply.  Note  that  we  already  used  the  term  demand  concerning  power  systems  and 
denoted  by  the  demand  appearing  in  the  area  that  has  to  be  satisfied  primarily  by 

system  j and  by  other  systems  only  in  the  case  if  system  j has  no  pov^r  enough  to  do 
that.  In  what  follows  will  be  termed  area  demand  (of  area  j)  or  simply  demand.  In 

our  planning  model  the  network  demand  at  system  j will  be  given  by  and  this  can 

be  positive,  zero  or  negative  as  well. 

If  S,T  are  subsets  of  N,  then  we  define 


q(S)  = I q(a)  , 
acS 

i(S,T)  = I i{a,by  . 
aeS.bCT 

Clearly  q(.S)  is  an  additive  set  function  and  i(S,T)  is  additive  in  S and  T whenever 
the  other  set  is  fixed.  More  explicitly,  we  have 

q(Sj^  bSj)  = q(S^)  + q(S^)  if  n = ♦ , 


i(S^  U S^,  T)  = 

i{Sj^,T) 

+ i{S2,T) 

if 

5^0  52  = 1). 

i(S,  Tj^  u Tj)  ” 

i(S,T^) 

+ i(S,T2) 

if 

Tf  n T2  = <> 
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These  definitions  iniply  for  every  S and  T 

i(S,S)  » 0 , 

i(S,T)  < y(S,T) 

The  network  demand  q is  called  feasible  if  there  exists  a flow  i such  that 
i(N,a)  > q(a)  for  every  a f N 

This  implies  that 

i(N,S)  > q(S)  for  every  S N 

Theorem  3.1  (Feasibility  theorem).  The  network  demand  q is  feasible  if  and  only  if  for 
every  S c N we  have 

(3.1)  q(S)  < y(S,S) 

where  S » N - S. 

In  the  above  presentation  it  was  more  convenient  to  write  i(a,b) , y(a,b)  and  q(a) 

instead  of  subscripted  variables.  Now  we  return  to  our  problem  of  interconnected  power 

systems  and  reestaJalish  the  notations  i , and  y The  notation  q(a)  or  q will  not 

ab  ab  ‘a 

be  used.  The  network  demands  are  5,  - x ,...,£  - x ; these  symbols  will  be  used  in  the 

XI  N N 

sequel. 

The  feasibility  theorem  provides  us  with  a system  of  linear  inequalities  expressing 
the  event  A(x,y,£).  The  inequalities  (3.1)  qive  necessary  and  sufficient  condition  that 
the  whole  pool  can  be  supplied  by  power.  The  y function  give  the  tie  line  capacities. 

Let  us  consider  some  examples.  The  first  is  a network  consisting  of  two  nodes  (see 
Fig.  1)  . Assume  that  y^^^  “ 72^  3et  them  denote  simply  by  y. 


Figure  1 . 


Then  q(l)  ” q(2)  = " ’‘2 


and  (3.1)  gives 


(S  is  empty)  , 


(3.3)  + i > Cj  . 

-y  < i < y • 

The  inequalities  (3.2)  give  the  solvability  condition  - with  respect  to  i - of  the  system 

(3.3) .  In  other  words,  (3.2)  is  the  projection  of  the  convex  polyhedron  (3.3)  given  in 

6 5 

R onto  the  R space  of  the  variables  x^,  x^,  y,  Cj* 

Consider  now  the  networ)i  of  three  nodes  illustrated  in  Fig.  2. 


1 


2 3 

Figure  2. 


Assume  that  the  tie  lino  capacities  are  symmetrical  and  Introduce  the  notations 
y . y ^ » y^j^j  yj  “ y23  ” ''32'  ''l  ’ ^^13  " '^31'  relations  (3.1)  give  the  following 

system  of  linear  inequalities 


■nie  convex  polyhedron  (3.4)  is  aqain  a projection  of  the  convex  polyhedron  (3.5). 

The  third  example  will  be  the  networ)c  illustrated  in  Fig.  3.  The  tie  line  capacities 


2 3 4 


Figure  3. 


i 
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are  assumed  to  be  symmetrical  and  will  be  denoted  by  y^  = ^12'  ^3  “ ^13'  '^4  ~ ^14' 
relations  (3.1)  are  specialized  as  follows.- 

^1  • ’‘l  ^ ^2  • ’'2  S ■ ’‘l  * ^4  - ^ ° ' 


Cl  - 


= ^2  ^ ^3  ^ ^4  ' 


^2  - ^^2 


63  - Xj 


«4  - ^ = 


?!  - *1  + S2  - *2 


Vj  ^ ''4  ' 


^1  - ’‘l 


+ £3  - X3 


= ^2  + ^4  ' 


5i  - x^ 


* «4  - *4  i ^2  ^ ^3  ' 


^2  ' *2  S ■ *3 


= ^2  ^ ^3  ' 


^2  - *2 


^ «4  - ^^4  = ^2  ^ ^4  ' 


S ■ *3  ^ ^4  - ^4  i ^3  * ^4  ' 


«2  - ’‘2  ^ ^3  - *3  * ^4  - ^ = ^2  * ^3  ^ ^4  ’ 


Ci  - x^  t - X2  f 53  - X3 


^1  - ’'l  * ^2  - ’'2 


^ ^4  - ^ ^ ^3  ' 


?!  - x^ 


t C3  - X3  + C4  - X4  J y2 


The  number  of  relations  contained  in  (3.1)  equals  2”-!  if  we  disregard  the  trivial 
relation  obtained  in  the  case  S ” N.  From  (3.1)  we  can  immediately  conclude 

Theorem  3.2.  The  linear  inequalities  determining  A(x,y,C)  consist  of  inequalities 


of  the  form: 


Partial  sum  of  the  ^j'®  5 


Partial  sum  of  the  x^'s  + Partial  sum  of  the  v -s  . 
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1 


Theorem  1 also  tells  us  what  partial  sums  have  to  be  combined.  We  can  write  them  auto- 
matically if  all  pairs  are  assumed  to  be  connected  with  eventually  zero  capacity  tie  lines. 
The  form  of  the  inequalities  emphasized  in  Theorem  2 will  be  important  when  proving  the  con- 
vex programming  property  of  our  optimization  problem  in  case  of  a multivariate  ganrna  distri- 


bution. 


4.  The  Use  of  Multivariate  Probability  Distributions 

A probability  distribution  defined  on  the  measurable  subsets  of  is  said  to  be 

logarithmically  concave  (logconcave)  if  for  every  pair  of  convex  sets  A,  B c r^  and 
0 < X < 1 we  have  the  inequality 

(4.1)  P(XA  + (l-X)B)  > tP(A)l^[P{B)]^‘^  . 

N 

If  the  rcindom  vector  C £ R has  logconcave  distribution,  then  the  same  holds  for  all  linear 
transforms  of  C [8].  Most  important  in  this  respect  are  the  following  theorems. 

Theorem  4.1.  If  C has  a continuous  probability  distribution  and  its  probadsility 
density  function  is  a logconcave  point  function  (i.e.  if  f denotes  the  density,  then  for 
every  u,  v e R*^  and  0 < X < 1 we  have  f(Xu  + (l-X)v)  > If(u)]^(f(v)  ^ ^)  , then  the 
probability  distribution  of  5 is  logconcave  i.e.  (4.1)  holds  where 

P(C)  = P(C  £ C) 

N 

for  every  measurable  C c r . 

Theorem  4.2.  If  gj^(z,w)  , . . . ,g^(z,w)  are  convex  functions  of  all  vari£d)les  in  the 
entire  space  and  5 is  a random  vector  having  (the  same  number  of  components  as  w and) 
a logconcave  probedillity  distribution,  then  the  function 

(4.2)  P(g^(z,e)  < 0,...,gj.(z,C)  < 0) 

is  a logconcave  point  function  of  the  variable  z. 

Theorem  4.1  was  proved  in  [7].  Theorem  4.2  easily  follows  from  Theorem  4.1.  For  a 
proof  see  [7] . Further  theorems  of  similar  type  are  proved  in  [2] . 

The  most  immediate  examples  for  logconcave  distributions  are  the  following: 
a.  Uniform  distribution  in  a convex  set  D c R^;  the 
probability  density  equals 


0 if  z ^ D , 


if  z £ D , 


where  D is  the  Lebesgue  measure  of  D assumed  to  be  finite: 
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b.  Normal  distribution;  the  probability  density  equals 


f(z)  = 


■y(z-in)  C ^(z-ra) 

, z £ r" 


where  C is  the  covariance  matrix  assumed  to  be  nonsingular  and  m is  the  expectation 
vector . 

We  shall  use  a third  multivariate  distribution  that  is  not  necessarily  logconcave.  The 
joint  probability  distribution  function  will  be  logconcave,  however.  This  is  the  multi- 
gamma  distribution  introduced  in  [10] . 

Let  ganma  distributed  random  variables  with  parameters 

^n'  ^0'  i = and  > 0,  j = 1, . . . ,N.  Thus  if  > 0,  then 

has  the  density 


^.-1  d,  -X.s 
X,3  z ^e  ^ 


z > 0 


If  d^  = 0,  then  by  definition  = 0)  = 1.  A gaista  distribution  is  said  to  be  stamdard 

if  the  X parameter  equals  1.  The  random  variables  are  clearly 

standardized.  Assume  that  every  is  a partial  sum  of  a collection  of  independent, 

standard  gamma  distributed  variables  nj^,...,n^,  where  r = 2”-!,  thus 

(4.6)  C - An  , 

where  A is  an  N x 2^-1  matrix  with  0,  1 entries,  all  columns  of  which  are  different 
from  each  other  and  from  the  zero  vector.  A representation  technique  producing  the 
(approximate)  (4.6)  form  for  5 is  given  in  (10). 

If  a standard  gamma  distributed  random  variable  n has  d parameter  > 1,  then  by 
Theorem  4.1,  n has  a logconcave  distribution.  If  d < 1,  then  write 


and  observe  t)iat 
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a.)  n is  a convex  function  of  n 


b.)  n has  a logconcave  density 


The  first  assertion  is  trivial.  The  second  assertion  can  be  proved  as  follows.  Let  z > 0: 

d 

then  the  probability  density  function  of  n at  the  point  z equals 

!/-> 

S’  d ^ 

. o / n < •*  \ = 

dz 


(4.8) 


1 

d S d S d 7 1 

35-  P(n  < ,)  = - P(n  < z ) = 


^i»-l  -t,^  1 

(tJ)  t e dt  p,d+i) 


1/^ 


The  function:  equal  to  zero  if  z ^ 0 and  equal  to  the  right  hand  side  expression  of  (4.8) 

d 

if  z > 0,  is  logconcave.  Thus  n has  a logconcave  distribution. 

Now  if  in  (4.6)  )ias  parameter  < 1,  then  we  represent  in  the  form  (4.7), 

otherwise  we  leave  n.  unchanged.  Then  express  C.  as  ■r5<  j=l,...,N  and  substitute 
1 j i 1 

the  thus  obtained  expressions  for  into  the  inequalities  determining  A(x,y,C).  Using 

Theorem  3.2  and  Theorem  4.2,  we  immediately  obtain 

Theorem  4.3.  tf  j = 1,...,N  and  C has  the  form  (4.6)  (where  n has 

^ j ^ 

independent  standard  gamma  components) , then 

P(A(x,y,C)) 

is  a logconcave  function  of  the  variables  contained  in  the  vectors  x and  y. 

As  an  exeui^le  consider  the  case  of  N » 2.  Assume  that  = 1,  - 1, 


'’2 


(4.9) 


^2  ' ’’>1 


n. 


where  « 2,  d^  « 1/2,  d^  * 3/2.  Substituting  and  intb  (3.2)  and  writing 

n,  » (•^)^  , 


we  obtain  the  inequalities 
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5.  Detailed  Formulation  of  the  Problem  for  N - 2 and  N - 3. 

In  case  of  N • 2 we  start  from  Relations  (3.3).  We  use  the  nonnegative  variables 
i,,  and  i but  omit  the  subscripts  for  the  sake  of  simplicity.  We  have  the  equality 


1-1  - 1 


The  case  that  at  the  same  time  i > 0,  i >0  will  be  automatically  discarded  by  the 
optimization  principle  provided  d^^  * 0 <^2i  * second  stage  problem  is  the 

following: 

minimize  ^21^~  ^1**1^  ^2^*2** 


subject  to 


xi  - i + i + j . 


*2  + i - i + ^2  - ^2  ' 


-y  < i - i ^ y 


i , i , Zj^,  Zj  J 0 


‘Rie  first  stage  problem  is  the  following  special  case  of  Problem  (3.1): 

minimize  Ic^(Xj^)  + ^^(x^)  + c(y)  + E(p)] 


subject  to 


tl  * ^2  ^ >'1  + >‘2 


p Ml  *y  \ip 


(i)  (u) 

X^'  < Xj  < ' 


Y*®’’  < y < y‘“’,  j - 1,  2 


where  y • *2'  ^1'  ^2^  (random)  optimal  solution  of  the  second  stage  pro- 

blem and  p is  a prescribed  (very  high)  probability.  The  12  subscripts  are  omitted  in 
c,  y and  Y<®’ , y‘“’ . 
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In  this  special  case  the  value  of  the  objective  function  of  the  second  stage  problem 
can  easily  be  expressed  in  terms  of  x^i  y.  Cj^«  £2"  cases  have  to  be  investigated. 

a. )  - Xj^  < 0,  < 0.  Then  vi  “ 0 . 

b. )  - Xj^  < 0,  ^2  “ *2  ^ optimal  solution  is 

i “ min(x^  - , y)  , i = 0,  = 0,  - x^  - i^  and  p = ‘^12^*  * *'2*^2*’ 

c. )  £1  ■ ’'j  £2  ” ’'2  - optimal  solution  is 

i = 0,  i = min(ej^  - x^,  x^  - V)  ' 22'°'  “^1  ° ^1  ” ’‘l  ' ^ “ ‘*21^'  *^l*’'l'' 

d. )  - Xj^  > 0,  ^2  ~ ^2  ^ optimal  solution  is 

” i “ 0,  Zj^  « £j^  - Xj^,  ^2  “ ^2  ~ ’'2  ^ “ ^1*^1*  * '^2**2*' 

If  £j^»  Cj  have  a continuous  joint  distribution  and  their  joint  probability  density 
function  is  f(Vj^,  v^)  . then  the  expectation  of  p can  be  written  in  the  following  meinner 

» ’'1 

E(P)  “ / / (dj^2  niin(Xj^  - v^^,  v^  - x^,  y)  + 


+ f2*''2  ’ ’'2  ” v^  - x^,  y)))f(Vj^,  v^JdVj^dv^  + 


/ ! (<^21  '"^"*''1  ■ ’'1'  ^2  ~ ''2'  * 


+ " ’‘l  ” x^  - v^,  y))]f(Vj^,  v^jdv^^dv^  + 


/ / - Xj^)  + ^2*'^2  ' *2’'^*''!'  ''2’‘^''i‘*'^2  ■ 


If  the  network  consists  of  three  nodes,  then  the  configuration  given  by  Figure  2 is 
the  most  general.  If  somewhere  there  is  no  interconnection  then  we  set  the  tie  line  capa- 
city equal  to  zero.  We  start  from  Relations  (3.5)  euid  use  the  representation 

^1  ” ^12  ■ ^12  ^13  ■ ^13  ' 

^2  " ^12  ■ ^12  * ^23  ■ ^23  ' 
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The  cases  > 0,  >0  will  be  automatically  discarded  in  the  problem  provided  the 

coefficients  of  these  variables  are  positive  in  the  objective  function.  The  second  stage 
problem  is  the  following: 


minimize Id^ji^j  + t d^3l*3  > d3^i-3  t d23i^3  + d33i;2  + 


r 


6.  Remarks  Concerning  the  Numerical  Solution. 


We  intend  to  return  to  the  numerical  solution  of  Problem  (2.1)  in  another  paper.  In 
the  special  case  when  E(u)  is  disregarded  in  the  first  stage  problem,  we  obtain  a proba- 
bilistic constrained  stochastic  programming  problem  similar  to  that  investigated  in  [11] . 
There  numerical  example  is  presented  using  the  normal  and  the  multigamma  distributions  and 
the  nonlinear  programming  method  of  supporting  hyperplanes  proposed  by  Veinott  [13] . Of 
course  other  nonlinear  programming  methods  can  be  used  too. 

When  computing  function  values  in  the  course  of  the  optimization  process,  the  LOLP  is 
obtained  in  case  of  the  assumed  capacities.  Thus,  as  a by-product  of  our  theory,  we  obtain  a 
method  for  finding  the  LOLP.  In  general  the  LOLP  are  obtained  by  simulation,  using  the  in- 
equalities defining  A(x,y,C).  In  some  special  cases  it  may  be  pxjssible  to  express  the 
LOLP  more  explicitly.  If  e.g.  C has  a discrete  distribution  of  some  simple  type,  then  we 
may  be  able  to  add  the  probabilities  belonging  to  those  C values  which  satisfy  the  inequal- 
ities A(x,y,C)-  This  sum  equals  1 - LOLP. 

Conclusion 

We  have  expressed  the  condition  that  interconnected  power  systems  can  assist  each  other 
to  meet  total  demand  in  the  pxjol  in  terms  of  lineeir  inequalities  containing  the  area  demands, 
area  generating  capacities  and  tie  line  connection  capacities  as  variables.  This  provides 
us  with  a tool  for  the  Monte  Carlo  simulation  of  the  system  reliability  = 1 - LOLP.  More- 
over, imposing  a lower  bound  on  the  system  reliability  and  minimizing  total  investment  cost 
we  are  able  to  plan  power  generating  and  tie  line  connection  capacities  of  interconnected 
power  systems.  A model  is  formulated  for  the  Sctme  problem  where  the  operating  and  outage 
costs  are  also  included.  In  this  case  the  niomerical  solution  is  still  possible  if  the  num- 
ber of  interconnected  systems  is  small.  The  method  of  solution  will  be  presented  in  a 
further  paper. 
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